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ABSTRACT 

We develop equations and obtain solutions for the structure and evolution of a pro- 
todisc region that is initially formed with no radial motion and super-Keplerian ro- 
tation speed when wind material from a hot rotating star is channelled towards its 
equatorial plane by a dipole-type magnetic field. Its temperature is around 10''K be- 
cause of shock heating and the inflow of wind material causes its equatorial density 
to increase with time. The centrifugal force and thermal pressure increase relative to 
the magnetic force and material escapes at its outer edge. The protodisc region of a 
uniformly rotating star has almost uniform rotation and will shrink radially unless 
some instability intervenes. In a star with angular velocity increasing along its surface 
towards the equator, the angular velocity of the protodisc region decreases radially 
outwards and magnetorotational instability (MRI) can occur within a few hours or 
days. Viscosity resulting from MRI will readjust the angular velocity distribution of 
the protodisc material and may assist in the formation of a quasi-steady disc. Thus, 
the centrifugal breakout found in numerical simulations for uniformly rotating stars 
does not imply that quasi-steady discs with slow outflow cannot form around magnetic 
rotator stars with solar-type differential rotation. 

Key words: circumstellar matter - stars: early type - stars: emission-line. Be - stars: 
magnetic fields - stars: rotation - stars: winds, outflows 



1 INTRODUCTION 

We use tlie term protodisc to refer to the disc region that is initially formed when wind material from a star is channelled 
by dipole-type magnetic fields towards the equatorial plane. The material in a protodisc region has no radial motion and 
rotates with the magnetic field. It is bounded above and below by shock boundaries. We develop models for the structure and 
evolution of protodiscs around magnetic rotator stars. We study protodiscs as separate entities while recognising that they 
may be precursors of quasi-steady discs, such as those of Be stars. We do not construct models of quasi-steady discs in this 
paper but indicate how protodiscs relate to models of Be star discs. For exa mple, a p rotodisc region could initially be present 
in the inner part of the quasi-steady disc model with slow outflow studied bv lBrown et al ( 2008,). Also, a more detailed stud: 
of protodiscs is imp ortant to verify how centrifugal outflow in protodiscs of rapidly rotating stars (e.g.. iMaheswaranI l200i 
ud Doula et a ] |2006l ) affects the formation of quasi-steady discs around magnetic rotator stars. In this paper we consider 
protodiscs around stellar models with both uniform and differential rotation. We study how their radial extents change with 
time and we determine the conditions under which MRI will occur to affect the angular velocity distribution and the radial 
outflow in the disc region. 

There has been a tremendous increase in observations of emission line B stars, including polarisation, interferometry and 
magnetic diagnostics, but the question of how discs with near Keplerian rotation speeds are formed around them is has not 
been fully answered. There c an be othe r mechanisms that are responsible for disc formation in different stars depending on 
their properties. For example, iLee et all (|l991i) have proposed a decretion disc model for Be stars with near maximal rotation 
speeds. There is no magnetic force in their model and angular momentum is transported by viscous forces. Be stars are 
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known to rotate rapidly (e.g., see lPorter fc Rivinius 2003 ) and the presence of significa nt magnetic fields in OB stars has been 
reported by several authors (e.g., 



Donati et al 



20061: lNeine3l200'7l : ISchnerr et ail2007l). The concep t of a Wind Compressed 
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Disc was initially developed bv lBiorkman fc Cassinellil ( 19931 ) . Subsequently. ICassinelli et all ( 2002h proposed a Magnetically 
Torqued Disk(MTD) model in which a sufficiently strong magnetic field channels a flow of wind material towards t he equatorial 
plane to for m a disc and provides the torque required by the disc to maintain quasi-Keplerian rotation speeds. MaheswaranI 
( 20031 . I2OO5I ) developed the Magnetic Rotator Wind Disk(MRWD) model in which Kepleria n discs may be form ed with stellar 
magnetic fields that do not have to be as strong as those required in the MTD model. Brown et a] 1 2004 ) presented an 
amended version of the MTD model to include stellar gravity darkening and lBrown fc Cassinellil (12005) suggested tha t this 
should be further modified to include the effect of gravity arising from finite thickness of the disc. Brown et al ( 20081 ) have 
studied the properties of discs with slow radial outflow and address criticisms of the MTD model. 

In the MRWD model for rotating stars with winds, when magnetic field lines that emerge from the stellar surface loop 
back to the star, a transient disc that has no meridional motion and corotates with the star is initially formed. After the 
local density in the disc has increased to a value where the centrifugal force excee ds the sum of gravi ty and the magnetic 
force, material diffuses radially outwards while stretching out the magnetic field lines. iMaheswaranI (|2003| ) suggested that such 
material can diffuse into quasi-steady discs and that the presence of viscosity durin g this stage would enable th e formation of 
Keplerian discs of larger radial extent. On the other hand, numerical simulations bv lud Doula et al (|2006l . l2008h for the flow of 
wind material along dipole magnetic field lines display centrifugal breakout in the form of transient bursts and they concluded 
that MTD type m odels were no t appr opriate to explain disc formation. However, there are problems with their simulations: (a) 
As pointed out by Brown et aJ ( 2008l l they did not follow changes for a sufficiently long duration of time, (b) The simulations 
are essentially two dimensional, (c) Numerical MHD simulations do not always provide accurate results because discretization 
of the MHD induction equation leads to the introduction of artificial magnetic diffusivity that is several orders of magnitude 
larger than the real magnetic diffusivity in discs. This can yield spurious results for magnetic reconnection that may not be 
realistic, (d) Their models fail to consider differentially rotating stars and, therefore, do not obtain differentially rotating disc 
regions that are subject to the onset of MRI. The occurrence of MRI can significantly affect the rotation speed in a disc 
region. We carry out an analytical study and obtain time-dependent solutions for uniform and differentially rotating stars. 
We find that when the stellar angular velocity increases along the surface towards the equator there is a broad range of stellar 
parameters for which MRI occurs in a sufficiently large part of the protodisc region. 

We consider a stellar model in which the magnetic field is dipole-type over some region of the stellar surface, so that field 
lines leaving the stellar surface pass through the equatorial plane and return to the star. Wind material flowing along such 
lines is channelled into a dense equatorial protodisc region. Since the flow of supersonic wind material into this region from 
above and below are in opposite directions, this material will cross shock surfaces that form the upper and lower boundaries 
of the disc region. We refer to the region between the sonic surface of the wind and the shock boundary of the disc as the wind 
zone. Thus, we need not be concerned about gas pressure effects on the wind flow, which is treated in an empirical fashion in 
that the velocity is chosen to obey a beta law and the density is assumed to follow from a global mass conservation relation. 
In section [21 we set up the relevant properties of axially symmetric wind zones that enable us to determine values of physical 
quantities at points on the wind side of the shock boundary from known values at the stellar surface. We develop formalism 
to account for the presence of differential rotation along the region of the stellar surface from which wind material flows into 
the disc region. The material in a protodisc region has no radial motion and rotates with magnetic field lines because the 
centrifugal and pressure gradient forces that are directed outwards are balanced by gravity and magnetic force. In section [3] 
we describe the general concept of a protodisc, which is a more generalised version of the pre-Keplerian disc developed in 
the MRWD model. Also, we describe the physical processes that are involved in their formation and evolution. In section 
[H we develop MHD generalisations of the Rankine-Hugoniot jump conditions that are appropriate for shock boundaries of 
protodiscs so that we can compute changes in values of quantities across shock surfaces. We take into account the vertical 
motion of shock surfaces caused by continuous filling up of a protodisc region by wind material. 

In sectionO we develop the equations that govern the structure of the models of protodiscs that we consider. We explicitly 
include the effect of the component of stellar gravity perpendicular to discs and, because of shock heating of the disc region, 
we include the vertical and radial components of the thermal pressure gradient in the respective momentum equations for 
protodiscs. Even though the heights of protodiscs may be relatively small, the presence of a vertical component of stellar 
gravity leads to a vertical gradient of pressure and density. The equatorial density of material in the protodisc increases with 
time as the disc fills up. This has the effect of increasing the centrifugal force in relation to the magnetic force as well as 
increasing the thermal pressure in relation to magnetic pressure. The first of these effects tends to cause material at the outer 
edge of a protodisc to move radially outwards whereas the latter leads to the onset of MRI within the protodisc. To determine 
the relative importance of these competing effects, in section [6] we derive criteria that influence the evolution of protodiscs. 
Also, we discuss the conditions under which the shrinking of a protodisc caused by centrifugal outflow is halted by the onset 
of MRI. This leads to viscous readjustment of the rotation speed of the disc region. In section [T] we apply our equations and 
jump conditions to protodiscs of different stellar models and present the results of our computations. In section [S] we discuss 
important factors that affect the evolution of protodisc material after the onset of MRI and refer to some relevant quasi-steady 
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models. Our conclusions regarding the formation and evolution of protodiscs around rotating magnetic stars with winds are 
summarised in section [9] An important objective of this paper is to relate the properties of protodisc regions to values of 
quantities at the stellar surface. Therefore, we include all the equations that are necessary to solve for the different physical 
quantities in the wind zone, at the shock surface and in the protodisc region. 



2 AXIALLY SYMMETRIC WIND ZONE 

In this section we present the equations that are used to determine the values of different physical quantities at points on the 
wind side of the shock boundary in terms of their values at the stellar surface. We assume that the wind zone is in a steady 
state and is isothermal. The mass, luminosity, radius and effective temperature of the star are AI, L, R and T^ff respectively. 
We use an inertial frame of reference with spherical coordinates (r, 9, (j)) in the wind zone and cylindrical coordinates {w, (j>, z) 
in the disc region. We assume that the star, wind and magnetic field are symmetrical about the axis of rotation. Let the 
angular velocity at any point P on the stellar surface be uj, which may vary with latitude along the stellar surface. We show 
in section |3] that stellar models with appropriate differential rotation will have protodiscs that satisfy a condition that is 
necessary for the occurrence of MRI. All the equations that we set up for the wind zone are valid even when uj varies along 
the stellar surface provided that the model is axially symmetric. We define an associated rotation rate parameter 5* for the 
point P by the equation uj = [GM / R^]^^'^ 5*. If the star has uniform rotation, then to is constant for the star and <S* gives the 
ratio of the equatorial rotation speed of the star to the Keplerian rotation speed. Otherwise, Si, will vary with lu along the 
stellar surface. Let p and B denote density and magnetic field, respectively. We use v and u, respectively, to denote velocities 
in the wind zone and the disc region. Subscript m denotes values of meridional components. Then, Vm ~ (wr +Vg^^^^ and 

Bm ~ (Sr + Bg^^^^. Let Q be a point of impact of wind material on the shock boundary. Subscripts W and D, respectively, 
denote values of quantities on the wind side and the disc side at the point Q. Let the meridional magnetic field line that 
emerges from the point P on the star pass thro ugh the point Q on the shock surface. 

Equations for the wind zone were set up in iMaheswaranI (|2003l ) using the steady state equations given bv lMestell (|l968l l. 



Initial values of different quantities in the wind zone are specified at the sonic surface and are denoted using subscript o. The 
MHD induction gives curl(?; X B) = 0, which implies that the meridional vectors Vm and Bm must be parallel to each other. 
Hence, the streamlines of flow in a meridional plane in the wind zone coincide with magnetic field lines. This, together with 
the mass and magnetic flux conservation conditions div(p?;) = and div 5 = require that 

P ^m Po '^m,o 



Bm B 

n 



= C, (1) 



where ^ is constant along meridional streamlines and represents the ratio of mass flux to magnetic flux. At the sonic surface, 
Vm,o will be equal to the isothermal sound speed. In general, ro and the density po at the sonic surface will be functions of 9. 
In our models, we take the sonic surface to be approximately spherical so that ro is constant. If the stellar mass loss rate M, 
is known only as an average value for the star rather than as a function of 0, we can write 

Anrlpo Vm,o = M , (2) 

where po is constant and represents the average density over the sonic surface. The azimuthal components of wind velocity 
and magnetic field satisfy the condition 

— LOUJ + VmB^i/Bm ■ (3) 

Equations ^ and Q allow us to visualise the magnetic field line that is anchored at the point P in the star to be rotating 
in the wind zone with the same angular velocity uj as P. The azimuthal field B^ is negative in the northern hemisphere of 
the wind zone so that field lines curve in the direction opposite to that of rotation. Since wind material flows along field lines, 
the rotation speed of the wind material is less than the corotation speed as verified in equation (|3]). In our stellar models, 
the material flowing into the disc region encounter the shock boundary before it could reach the Alfven point. The azimuthal 
torque equation for the wind can be written in the form 

Bm.o B^ \ 

'"<t> ~ = -C, (4) 

'i-KpoVm,a J 

where C is constant along a streamline and may be interpreted as total angular momentum per unit mass. Each field line in 
the wind zone is associated with a characteristic density pA and a characteristic distance zua from the rotation axis. These 
are given by 

" ~B^ " 4^r4 B^ ^ ' 

and 
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zua ^ — = [ 1 ~ • (.o; 

LO \ LJ / y iivpo Vm.o V^^o J 

In general, the values of pA and vja will vary with the latitude of the point P on the stellar surface from which the corresponding 
field line emerges. Note that a streamline that does not enter the disc region but travels with the wind to a large distance 
will have an Alfven point where p = pA and ro = zua ■ However, for streamlines that enter the shock boundary of a protodisc 
region, there is no point in the wind zone where p and -co are equal to these Alfven values but we define them as streamline 
constants because we can use them to compute the rotation speed and azimuthal component of the magnetic field at any 
point in the wind zone. We solve equations (|3]) and Q to obtain the values of the ^-components of the wind velocity and 
magnetic field. Then, using equations ([5]) and ((6)1 we obtain 

, ^ , . [l - {pat^a) I (P^^)] (1 - PaIpo) \ 
v^,w = v^,o — < -iq , 2 ^ „2M TTTT } (V 




[1 - {pAZ^\) I {p^vjI)\ (1 - paIp) 



and 



_ {l-^\l^^)[l-pAlp.) \ 




We express the meridional magnetic field in the wind zone in the form for dipole-type fields developed in [Maheswaran 



I 2003h . We have Br = 2C (ro/r)* cos Q and Bg = {b - 2)C {ra/rf sinO, where 6 is a constant such that 2 < 6 ^ 3. We should 



point out that 6 = 3 represents a strict dipole field and that the the meridional component of the magnetic force will be zero 
in this case. However, for 2 < fe < 3 there is a nonzero meridional component of the magnetic force. As was pointed in the 
earlier paper, even if 6 = 3 at the stellar surface its value will be less than 3 in the wind zone because the field lines will be 
stretched out to some extent by the wind material. We obtain the meridional component of the magnetic field in the wind 
zone as 



Bm ~ C 



(r^)' [4-6(4 -b)sin^e]'/' . (9) 

In the case of disc regions with relatively small height, we have sin 6 « 1 at points along the shock boundary. Then, equation 
((9| gives an approximation for Bm,w in the form 

B^,w^{h~2)c{^y . (10) 

To obtain the meridional speed Vm,w of the wind material as it enters the shock boundary, we use the beta velocity law for 
line driven winds given by 

Vm,W ^Voc[l-jj , (11) 

where /3yei is a constant that is in the range 0.5 < /3„ei < 2 (see e.g., Lamers fc CassineUi 19991 '). Also, T4o is a terminal speed 



that is constant for each streamline but may be different for different streamlines. In our computations we take Voo to be the 
same for all streamlines and equal to the velocity at infinity of the wind. The density pw of the wind at the shock surface is 
obtained in terms of the density at the initial surface given by equation ([2]) by using the mass and magnetic flux continuity 
equation ([T} together with the velocity law pip . In stellar models where the ratio F of the stellar luminosity to the Eddington 
luminosity becomes important, we replace M by M(l — F) in the appropriate equations. 

We introduce the variable x = ru/R to denote the non-dimensional distance x of the point X on the equatorial plane, whose 
distance from the centre of the star is ru. When a subscript is added to X we add the same subscript to the corresponding x. 
The wind material that is channelled along flux tubes of a magnetic field crosses a shock surface and moves into an equatorial 
region between the stellar equator and a point Xum, which is the outermost point on the equatorial plane to which the 
magnetic field is able to channel the fiow of wind material. Let Xkep be the innermost point on the equatorial plane where 
local Keplerian rotation speed equals the corotation speed. Then, x^ep = S^^^^, where Sj, is the rotation rate of the star at the 



point Pkep on the stellar surface that is linked to Xk^p by a magnetic field line (ICassinelli et all 120021 '). The conditions under 
which the stellar magnetic field will be abl e to channel the win d into a disc region and relevant parameters are discussed in 
detail in sections 3, 4 and 10 through 13 in iMaheswaranI (|2003l '). In particular, the value of the meridional Alfven paramet' 



er 



Am ~ Bml{\/^T^p Vm) must be larger than unity at all points along a streamline from the star to the shock surface. Hence, 
the strength of the meridional magnetic field at the initial surface must be larger than a minimum value given by Bm.o,min 
for the formation of protodiscs. Values of Bm,o,min are approximately IG for B5 or B9 stars, lOG for B2 stars, 60G to 80G for 
BO stars and of the order of lOOG for O type stars. Our interest is in stars with magnetic fields that are only slightly stronger 
than these minimum values so that they are able to channel wind fiow into a protodisc region. In this paper we consider 
models of O and B type stars shown in table [T]with rapid rotation such that 5* > 0.5. 
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Figure 1. Schematic figure of a meridional cross-section of a protodisc region of the star. Dashed fines represent the coincident streamfines 
and magnetic field lines. 



3 PROTODISC REGIONS 

Here, we describe the physical properties of protodiscs and the processes involved in their formation and evolution. A magnetic 
field of sufficient strength in a hot rotating star channels the flow of wind material from the star into an equatorial region. 
Furthermore, in some part of this region, the field lines initially loop back to the star and the magnetic force is strong enough 
to provide support against centrifugal force so that the material accumulates to form a protodisc with no meridional motion. 
It is bounded above and below by shock surfaces as a result of the fiow of wind material towards the equatorial plane from 
opposite hemispheres. Figure [1] gives a schematic meridional cross-section of a protodisc region in an axially symmetric star. 
Let 5 denote the angle that the tangent line to the shock surface at the point Q makes with the equatorial plane. We use 
subscript n to denote the normal component of a vector at the shock boundary in the direction from the wind zone into the 
disc region. Let h be the height of the point Q above the equatorial plane. We assume that the material in the disc region 
satisfies the ideal gas law p = TZpT/fi, where TZ is the gas constant and fi is the mean molecular mass. For convenience, we 
use the substitution p = pa^, where = TZT/p. Note that the adiabatic sound speed is a'y^^^, where 7 is the usual ratio of 
specific heats. 

Let Xinner and X^dge showu iu Figure [1] be the innermost and outermost points of the region in the equatorial plane 
where the magnetic force is strong enough to prevent radial motion against the effects of the centrifugal, gravitational and 
thermal pressure forces. The point Xi„„er will be located very close to the point Xkep and its position will be infiuenced by 
the radial components of the magnetic and thermal pressure force in addition to gravity and centrifugal force. A protodisc will 
form in the region between the points Xinner and Xedge- Wind material arriving between the stellar equator and Xi„ner will 
tend to fall back towards the stellar surface because the sum of gravity and magnetic force dominates the sum of centrifugal 
and pressure forces. Material arriving between X^dge and Xum will fiow outwards because the centrifugal force dominates in 
this region. Meridional magnetic field lines emerging from the points Pinner and Pum on the stellar surface are connected to 
the points Xinner and Xum, respectively, on the equatorial plane. 

The protodisc region has super-Keplerian rotation speed. As this region fills up, the density at points along the equatorial 
plane increases and the centrifugal force increases relative to the sum of gravity and magnetic force. While the location 
of Xinner rcmaius fixed, centrifugal outfiow causes the outer edge of a protodisc X^dge to move towards Xinner as time t 
increases, provided that factors such as MRI do not intervene. At the same time, the continuing impact of wind material 
on the shock boundary causes the protodisc to heat up. The thermal pressure can rise to values where it exceeds the local 
magnetic pressure and the ratio fimag ~ 9,11 pa? /B'^ can become larger than unity. Let the angular velocity of the disc material 
at a point X on the equatorial plane be Q. so that u^,x = O.'co. There will be an onset of MRI when the two conditions 

do. 
dvj 



Prnag > 1 and — < (12) 



are satisfied (jBalbus fc Hawlevlll99ll ). Then, viscosity will become an important factor in the rotation of the disc region and 
the magnetic field will be affected by turbulence resulting from MRI. 

In general, the angular velocity u of the region of the stellar surface from which material flows to a protodisc region will 
vary with latitude. From section O we know that lo is constant along a meridional magnetic field line in the wind zone. In our 
model of a protodisc, each magnetic field line from the stellar surface that enters the protodisc region is continuous across 
the protodisc and returns to the stellar surface. In the case of stars that have differential rotation in which u) increases along 
the stellar surface from the point Pum to the point Pinner, if the shock boundary of the protodisc is stationary, the protodisc 
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rotation theorem established in iMaheswaranI (|2003h shows that there wiU be a corresponding decrease in the angular velocity 
Q. of the protodisc region from Xinner to X^dge- That is, the second of the two necessary conditions (|12p . dQ./dzo < 0, for 
the onset of MRI in protodiscs is satisfied when the angular velocity ui increases along the stellar surface towards the stellar 
equator and the shock boundary of the protodisc is stationary. In section |6l we consider the case where the shock boundary 
of the protodisc is not stationary but moves so that its height changes with time. Here too, results of numerical computations 
presented in section [7] show that for appropriate models of differential rotation of the star where lj increases towards the 
equator, the condition dil/dzzi < is satisfied in the protodisc. The first condition for the onset of MRI wil l be satisfied when 



the in crease in density and shock heating of the protodisc region cause Pmag to become larger than unity. iBalbus fc Hawlev 



1 19911 ) have also shown that the rapidly growing modes increase exponentially on a time scale of order 4/(3 0) seconds. The 
rotation time scale for the protodisc region is trot = 2tt/Q seconds. Thus, MRI will occur on a time scale that is even shorter 
than the rotation time scale when the necessary conditions for the onset of MRI are satisfied. Although we focus on the onset 
of MRI, it is possible that other types of MHD instability could occur in protodiscs. Our interest is in identifying at least 
one mechanism that will prevent a protodisc from shrinking completely and redistribute its angular momentum so that a 
quasi-steady disc may be formed. 

We consider the protodisc phase of a disc to end at the time that MRI occurs in that region. In an axially symmetric 
protodisc, the meridional component of the velocity is negligible and we have u = {0,u^,0) in cylindrical coordinates. The 
magnetic field is given by i? = (Bin, B^, Bz). When h/zu is relatively small, the slope tan 5 of the shock surface is also small 
and the normal component wind velocity v„,w is a function of only along the shock boundary. Then, because the ram 
pressure from inflowing wind material on the shock boundary remains steady, the temperature To of the protodisc at a point 
Q of the shock will also be function of zu only and does not change as h increases. Thus, if the radiative cooling time of the 
protodisc is much longer than the flow time of wind material into the shock surface, the temperature T in the protodisc will 
have the same value along a vertical line because there is no meridional flow of material in a protodisc region. That is, we can 
take T to be independent of z so that T = T{zu) = Toijii) and da/dz = 0. In models where h/u} becomes a large fraction 
before the onset of MRI it will be necessary to solve for T as a function of zu and z. 



4 JUMP CONDITIONS ACROSS PROTODISC SHOCK BOUNDARIES 

In this section we set up the equations that are necessary to determine values of quantities at points on the disc side of 
the shock surface from the values on the wind side that we obtain in section [2l As material flows from the wind zone into 
a protodisc region, the height h of the shock surface forming the boundary of the protodisc increases. Let —Ua denote the 
normal component of the velocity of the point Q on the shock surface in the direction of n. To solve for the properties of a 
protodisc we need the values of po, Pd, Us, U0,d, Bn,D and B^^d at points on the disc side of the shock boundary. These are 
obtained by solving the six Rankine-Hugoniot jump conditions (|13|l through (|18|) . 
The continuity of mass across the shock boundary gives 

PW {Vn,W + Us) = PdUs ■ (13) 

Because the normal component of the magnetic field and the tangential component oi v y. B must be contiimous across the 
shock boundary, we have 

Bn.W = Bn,D (14) 

and 

{Vn,W + Us) B^^w — V^^w B„,W ~ UsB^^o — 'U;j>,DBn,D ■ (15) 

Continuity of the normal component of momentum gives the jump condition 
B^ B^ 

pw + pw {vn.w + Us) + —r^ = PD + pnUs + -^^ (16) 
and the equation for the continuity of the 0-component of momentum is 

, , rr ^ B„,wB,f, w B„^dB^^d /,_s 

pw (yn,w + Us)vct,,w -. ^ PdUsU4,,d z ■ (17) 

47r 47r 

The jump condition for energy across the shock surface is obtained by using the condition that the sum of the kinetic energy, 
internal energy, thermal energy and electromagnetic energy must be conserved. Then, 
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for the flow of energy per unit area normally across the shock surface. When both sides of this equation are expanded to 
remove the brackets, we obtain four terms on each side. The first term represents the flow of kinetic energy of the material 
and the second term is the sum of the flows of internal energy and thermal energy. The sum of the third and fourth terms 
gives the normal component of the Poynting vector that represents the flux of electromagnetic energy of the plasma. The 
Rankine-Hugoniot equations that we have set up above have general validity for protodiscs. However, in the case of the 
protodisc models that we consider here, the changes in values of and across the shock surface are not large. Thus, we 
flnd that the influx of mass, momentum and kinetic energy of the wind material go mainly into increasing the the density and 
temperature of the protodisc. 

Equations (|15|) and (|17p give the azimuthal components of velocity and magnetic fleld on the disc side of the shock surface 

as 

Vn.W B„.W Bs w „\ 

and 

where An,D = Bn.D / \/47i"PD is the normal component of the Alfven velocity at a point on the disc side of the shock boundary. 
In our protodisc models we have An^o larger than Us ■ This is a consequence of the condition that the stellar magnetic fleld 
is able to channel the flow of wind material into the disc region. Since -B^,^ < 0, equations (I19|l and (|20[1 show that the 
rotation speed and azimuthal component of the magnetic fleld increase across the shock boundary. In our protodiscs we flnd 
that u^^o is at most 2% larger than the corotation speed lu-uj. In protodiscs where the shock s urface is stati o nary, that is 
Us = 0, we flnd that u^^o = ujzu, which is consistent with the protodisc rotation law given in Maheswaran 1 20031 ). Also, 
because I-B^.ijl < |-Bc/),w|, the strength of B^ decreases slightly across the shock boundary. 



5 EQUATIONS FOR PROTODISC REGIONS 

Here, we set up equations for the structure of protodiscs and obtain analytic solutions. We do not consider a disc region to be 
a protodisc after the onset of MRI. Since material flows continuously into a protodisc region along magnetic flux tubes and 
does not flow outwards, its height h increases with time t. The accumulation of mass can be expressed in the form 

pdz = pwVn,wt- (21) 

Let X be a point on the equatorial plane at distance zo from the centre of the star and let Q be on the shock surface 
vertically above X. To obtain values of the magnetic fleld components at the shock boundary, we assume that during the 
protodisc phase h is relatively small compared with tj} and that the angle 5 is small. The computational results that we obtain 
for our models conflrm that this assumption is appropriate. Then, we have Stu.D = Bn.D sin 5 ~ and -Bz,d = — -Bn,D cos 5 « 
—Bm,w- In such protodisc regions, the meridional magnetic fleld, which is a continuation of the dipole-type fleld in the wind 
zone, is almost vertical and we may take Bti, ~ and B^ ~ —Bm,w- 

The general form of the azimuthal component of the momentum equation for an axially symmetric disc region is 

^ = -Mti7-7r-^ - Uz— ^ ^ + 2 75— ( a p c J , (22) 

ot ozu oz zu p pvj-' azu ^ ' 



where a is the viscosity coefficient developed by Shakura fc Sunvaev ( 19731 ). Cs is the local sound speed and 
1 



B^ ^ {ZUB4,) + -Bz 1^ {ZZB4,) 

ozu ^ ^' (jz 



(23) 



is the 0-component of the magnetic force. Before the onset of MRI, material in a protodisc region and the magnetic field 
rotate steadily with super-Keplerian speed and a = 0. Since u = (0, it^,0), we require solutions that satisfy the conditions 
du^/dt — 0, dB/dt — and divB = 0. Then, equation (|22p requires that dB^/dz = 0, which says that the azimuthal 
component of the magnetic field between the equatorial plane and the upper shock boundary is constant along a vertical 
line. Applying the boundary condition at the shock surface we have B^ = B^^o- Note that, in the part of the protodisc 
below the equatorial plane, B^ will be positive and have the same magnitude as in the upper part where it is negative. This 
implies a discontinuity in across the equatorial plane, which is physically plausible and is accounted for by the presence 
of a current sheet. The condition that dB/dt — together with the MHD induction equation imply that curl(u X B) — 0. 
When we substitute the solutions that we have already obtained for B^, -B^ and Bz we find that satisfies the condition 
du^/dz = 0. Therefore, as in the case of B^, the rotation speed is the same at all points along a vertical line in the 
protodisc and using the boundary condition at the shock surface, we obtain — u^^d- Thus, we have a physically valid 
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model of a protodisc that is described by the self-consistent solutions u = (0,u^,d,0) and B = {Q, B^^b, —Bm.w)- Also, we 
find that the vertical component of the magnetic force JT^ = in the protodisc region of our model and the radial component 
T^n.x ~ (6 — 1) -B^ I^-KVJ at the point X on the equatorial plane. Since u-^ is negligible in the protodisc region, the general 
form of the z-component of the momentum equation for such a region is 

dv pGM z 1 



dz r i-K 



0B^_8B^^^^^(dB, 



dz duj J \ d. 



(24) 



where p and p are functions of time t. This equation implies that the pressure gradient in the a-direction at any point is 
instantaneously balanced by the z-components of stellar gravity and magnetic force. In the case of our protodisc models in 
which z is much smaller than r we may take z/r'' ~ zj-uj'^ . Since p = pcP" and a is independent of z, equation (|24p gives 

dp _ 2z 
Yz ^ ~li 



2P- (25) 



where H = H{zu) — {2zu^a? /GM'^^^'^ is the density scale height in the protodisc. This equation implies that the density 
gradient of the protodisc in the vertical direction is influenced by the vertical component of gravity. The vertical component 
of the magnetic force is zero in our model and does not contribute to vertical equilibrium. We solve this equation to write the 
density of the protodisc at height z above X, in the form 

p = PDe^ , (26) 

where po is the density of the protodisc at the point Q of the shock boundary. Using equation (|21|l and the expression for p 
in equation (|26p we obtain an equation connecting the height h of the protodisc and time t in the form 



H 

t 



pWVn,W 



^.p,e erf^-j 



(27) 



where Bri{h/ H) denotes the error function value at h/H. 

Equations (|26|) and (|27p show that the density at any point X on the equatorial plane of the protodisc region increases 
with time. Then, unless MRI intervenes and terminates the protodisc phase, the density at X will increase to the value px,edge 
that is required for the centrifugal force to balance the sum of gravity, magnetic force and thermal pressure. That is, px,edge 
is obtained as the solution of the equation 



4 _ GM\ _ d (pa^) _ (6 



Px.dg. -5 - - - - 'S- - . (28) 



J dzu 47r zj 



We denote the time at which this occurs by tx,edge and let hx,edge be the height of the protodisc at X^dge at this time. The 
value of hx,edge IS obtained by putting p — px.edgc in equation (|26|) and solving for h. The time tx,edge at which this occurs 
is found by using equation (|27|l . In addition to giving us the values of the properties at the outer edge, the importance of 
these equations is that they give the rate at which the radius Xedge of the outer edge of the protodisc decreases. Thus, we are 
able to compute the radial extent of a protodisc at the time at which MRI occurs. Numerical values are shown in the tables 
presented in section [7] 

In our model of a protodisc, a magnetic field with a sufficiently strong meridional component channels the flow of wind 
material into an equatorial disc region. A much weaker azimuthal component of the magnetic field acts in concert with the 
meridional component to apply a torque on the wind material so that it enters the protodisc region with super-Keplerian 
angular velocity after crossing its shock boundary. The radial component of the magnetic force in the protodisc is initially able 
to assist gravity in maintaining radial equilibrium against the centrifugal force and thermal pressure. The protodisc region 
rotates steadily with no meridional motion when and B^ are constant along vertical lines. Consequently, the vertical and 
azimuthal components of the magnetic force are zero in the protodisc region and the density gradient in the vertical direction 
is influenced only by stellar gravity. Increase in equatorial density causes centrifugal outflow at the outer edge and causes the 
radial extent of the protodisc to decrease. It will continue to do so unless MRI or some other MHD instability occurs. We 
discuss this in section |6l 



6 ONSET OF MRI IN PROTODISCS 

As we pointed out in sections [3] and |4l continuous flow of wind material causes contiimous shock heating of the protodisc. 
The inflow time of the wind into the protodisc region is 

iinflow ~ ! TT" . (29) 

' V„ + Us 

The temperature of a protodisc will remain high if the shock heating is not offset by radiative cooling. The radiative cooling 
time for the disc material is given by tcooi = U/Lcooi ~ 3kT/{2NeA) s, where U = SkTNtotaifi ergs/cm"^ is the internal 
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Table 1. Properties of the different stellar models considered 



Star Spectral M R T^ff M Vca Pa Vm.o '",p,-k,eq,crit 

Name Type {Mq) (Rq) (lO'^K) {fQ-^Mgyr-i) (km s'^) (gm cm'^) (10''cms-i) (km 8"^) 



Model 06.5 


06.5 


29 


10 


4.0 


310 


2500 


1.38 X 10-12 


2.32 


661 


Model BO 


BO 


15 


6.6 


3.2 


27 


1300 


3.09 X 10~13 


2.08 


634 


Model B2 


B2 


8.3 


4.5 


2.3 


0.4 


840 


1.16 X 10"" 


1.76 


589 


Model B5 


B5 


4.5 


3.5 


1.5 


0.01 


580 


5.95 X 10-1^ 


1.42 


495 


Model B9 


B9 


2.6 


2.6 


1.0 


0.0013 


460 


1.72 X 10-1*^ 


1.16 


437 


a Ara 


B2 


9.8 


4.8 


2.3 


0.45 


1540 


1.15 X 10-1-* 


1.76 


625 


a Eri 


B3 


6.9 


4.1 


1.9 


0.042 


1330 


1.63 X 10-15 


1.62 


569 


a Col 


B7 


3.7 


3.4 


1.3 


0.003 


1250 


2.06 X 10-1^ 


1.32 


457 



Entries for the stellar properties are from lBiorkman &: Cassinellil lll993h for the theoretical models and from lCohen et a] (Il997l) for 
specific stars. Here, f 0,*,e9.c7 it is the critical Keplerian rotation speed at the equator. 



energy per unit volume and Lcooi = Ntota A erR s/(cm^s) is the cooling rate per unit volume, where A'^ denotes the 



number of particles per cubic centimetre. ICoxl l|2000f ) gives log A = —21.65 — 0.6 (log T — 5.5). From this we obtain the 



expression ergs cm /s for the cooling function that is appropriate for the postshock temperature ranges 

of our protodiscs. We assume that the disc material consists of Hydrogen and Helium with mass fractions Xh and 1 — Xh, 
respectively. Then the number density of electrons is given by A^e — (1 + Xh)pd / {2mH) per cubic centimetre and the cooling 
time is 

_ 6.72 X IQi^ k ruH Tj^^ _ 1.55 x IQ-^i T^*^ 

" {l+X„)pn " {1+Xh)pd ^ ' 

Our computations show that the radiative cooling rate in our models is slower than the rate of shock heating so that protodiscs 
continue to maintain their high temperatures. 

There are several values of the density and time scales that play critical roles in the evolution of protodiscs. The require- 
ment Pmag > 1 for the onset of MRI stipulated in condition (|12|l will be satisfied when the density px at a point X on the 
equatorial plane exceeds the critical value px,MRi given by 

P""''""' ^ W = ' ^^^^ 

because B^^x = 0. If Px.mri < Px.edge at a point X in a protodisc, the value of px will reach the value of the density required 
for MRI to set in at X before the outer edge of the protodisc reaches the point X. Let tx,MRi denote the time required for 
MRI to set in at the point X. 

Let the point P on the stellar surface be linked to the point X on the equatorial plane of a protodisc region by a meridional 
magnetic field line. Let lu be the angular velocity of the stellar surface at P and let be that of the protodisc at X. In section 
O we pointed out that in models of protodiscs where the shock surface is stationary, the condition that uj should increase 
towards the equator along the stellar surface will ensure that the second necessary condition, dQ,/dzo < 0, for the onset of 
MRI is satisfied. However, in our models of protodiscs, the shock surface moves so that its height increases with time. In this 
case, equations ([Sjl, (I14p and (|19p give 

n^uj ^-,wB^,wU^s (32) 

mB^,w{Al^o-Ul) 

where we have used the equations Q,-uj = = u^^d- We interpret th i s to b e a generalised form of the protodisc rotation law. 
Note that we recover the protodisc rotation theorem in iMaheswaranl \2mi ) when Us = 0. For stars with difi^erential rotation. 



we need the distribution of to along the stellar surface to compute values of fi. To determine whether will decrease with 
increasing in a protodisc for realistic values of oj, we use the angular velocity model for the Sun (e.g., see Solar Rotation 
in Allen's Astrophysical Quantities 2000), which has the form 



0=' 



1 - e 



Ab-2) 



(33) 



where u^q is the angular velocity of the star at the equator, denotes colatitude of the point P on the stellar surface, x is 
the non-dimensional distance of the point X from t he rotation axis an d e is a constant such that ^ e < 1. We have used 
the relation sin^ 5* = l/x^^~'^'' for dipole-type fields ( Maheswaran 2003h . When uj varies along the stellar surface, so does the 



rotation rate 5*. We use 5*, eg to denote the value of Si, corresponding to ujeq. In section [71 we present results of numerical 
computations for a variety of models. We find that the condition dQ,/dw < is satisfied in protodisc regions for a wide range 
of permissible values of e. 
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Table 2. Variation of angular velocity Q and values of X2 for protodiscs of differentially 
rotating stellar models with the parameters of a Ara, dipole 6 = 3, /3^ei = li Sir^eq = 0.7, 
Ldeq = 1.312 X 10~^ radians per second and e = 0.2 





= 57G, X2 = 2.31 




= 119G, X2 = 


2.37 




= 245G, X2 = 


2.45 


X 


n 


X 






X 


n 






radians/s 




radians/s 






radians / s 




1.31 


1.250 X 10-4 


1.31 


1.250 X 10" 


-4 


1.31 


1.250 X 10- 


-4 


1.61 


1.213 X 10-" 


1.61 


1.213 X 10" 


-4 


1.61 


1.213 X 10" 


-4 


2.01 


1.195 X 10-4 


2.01 


1.183 X 10- 


-4 


2.01 


1.181 X 10" 


-4 


2.31 


1.192 X 10-4 


2.31 


1.168 X 10" 


-4 


2.41 


1.159 X 10" 


-4 


2.41 


1.194 X 10-4 


2.41 


1.165 X 10" 


-4 


2.51 


1.156 X 10" 


-4 



7 RESULTS FOR PROTODISC MODELS 

We have applied the structure equations and shock jump conditions to study protodiscs around a variety of stellar models 
with properties listed in Table [l] We have performed computations for uniform and differential rotation with different stellar 
rotation rates and different surface magnetic field strengths. We have used equation (|32|l to compute values of the angular 
velocity Q of protodiscs around stellar models with solar-type differential rotation given by equation (I33p . for different values 
of e. In these models, the stellar angular velocity lj increases outwards along the surface towards the equator. We display 
numerical results for protodiscs around stellar models with uniform rotation and for solar-type differential rotation with 
e = 0.2. In section[2]we mentioned that b will be less than 3 in the wind zone. Here, we show results for b — 3 because they are 
not significantly different from those for models with values of b slightly less than 3. Table [2] shows the values of the angular 
velocity fl of protodiscs around differentially rotating stellar models with the parameters of a Ara, dipole b — 3, Pvei = 1, 
5*,e<j ~ 0.7 and ujeq = 1.312 x 10~* radians per second. These results show that dfl/dzu < for values of x between Xi„ner 
and X2- We find that in models with uj given by equation (|32[> and e > 0.1, the condition d^l/duj < for the onset of MRI 
is satisfied. In the case of a stellar model with uniform rotation where uu is constant along the surface, our numerical results 
show that the value of the angular velocity fl of the protodisc is constant to an accuracy of 0.1% or better for B-type stars 
and for O-type stars the variation in the value of Q is of the order of a few percent. 

In all our models, the outer edge Xedge of the protodisc initially coincides with the point Xi where the sum of the 
magnetic force and gravity is equal to the centrifugal force. The rotation speed at points along a protodisc is super-Keplerian 
before the onset of MRI. As material continues to flow from the wind zone into the protodisc region, the density along the 
equatorial plane increases and material at the outer edge starts to flow away radially so that X^dge moves towards the inner 
boundary Xin„er- In models with uniform rotation, the protodisc will continue to shrink radially unless some MHD instability 
intervenes. The rate at which X^dge approaches Xinner becomes slower as it gets closer to Xi„„er. In some of the cases the 
time taken for Xedge to approach Xinner can be many years. We use tahrink to denote the value of tedge when a protodisc 
is not affected by any instability and shrinks to a thin ring. In such a case, the radial extent of the protodisc will be small 
during most of this time and the protodisc will have the appearance of a ring. Our numerical results confirm that the values 
of h are small compared with zu during the protodisc phase in all the models that we consider. 

In a protodisc where the angular velocity uj of the central star increases outwards along the stellar surface towards the 
equator, we find that as Xedge moves inwards, it reaches a point X2 at the same time t2 at which MRI sets in at X2. At all 
points between Xinner and X2 we find that MRI sets in before the local density increases to the value required for centrifugal 
outfiow. We do not follow the evolution of disc material after the onset of MRI when viscosity becomes important. In some 
cases where MRI sets in after the protodisc has shrunk considerably, a quasi-steady ring may be formed. 

Table |3] shows the radial extents of protodiscs and the locations and times at which MRI occurs in our protodisc models. 
The values of Xi„ner, xum, xi and X2 and the time t2 are shown for protodiscs of different stellar models with uniform 5* 
equal to 0.5 and 0.7, and for differential rotation with 5*, eg ~ 0.7 and e = 0.2. In most cases, MRI sets in before the protodisc 
shrinks to a thin ring and occurs on a time scale of hours or days. In situations where the radial extents of protodiscs decrease 
to small values, we find that the time taken may vary from several days to a few years depending on the type of star, its field 
strength and its rotation rate. In all our models, the protodisc temperature To is of the order of 1 to 10 x lO^K and the wind 
inflow time tinfiom has values in the range from about 1 minute to about 40 minutes. We use equation (|30[) to derive the 
expression tcooi = 1-82 x 10~^*T^'®/p seconds for the radiative cooling time for a plasma with 70% Hydrogen and 30% Helium 
by mass. In all our models, we find that the values of teooi have values in a range from about 10 hours to several days. Thus, 
tcooi is many orders of magnitude larger than the inflow time tinfiom of wind material. Also, in all our protodisc models, the 
values of tcooi are much longer than tedge when x is large. Thus, the outer parts of protodiscs do not cool down before the 
material in those parts flow away from the protodisc. The speed Uahock of the shock boundary is of the order of 10 to 100 
kilometres per second and depends on the stellar model and the radial distance from the axis rotation. In all the protodiscs. 
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Table 3. Properties of protodiscs of theoretical stellar models listed in Table [T] with uniform rotation or differential 
rotation and aligned dipole-type magnetic fields 



Uniform Rotation with cS* = <S*,e(j Differential Rotation with e = 0.2 



Star 
Model 




(G) 




(G) 




Xl 


^shrink 

(days) 


(G) 




Xl 


^2 


t2 

(days) 


06.5 


0.50 


194 


1.78 


-26.81 


1 

1. 


.oy 


1.74 


0.03 


-24.47 


1 

1, 


.Do 


1.77 


1.76 


0.006 






432 


2.61 


-36.09 


1 

1. 




2.36 


0.22 


-31.64 


1 

1, 


.Do 


2.43 


2.42 


0.007 




0.70 


194 


1.78 


-37.54 


\ 


.27 


1.58 


0.18 


-34.25 


1 


.31 


1.62 


1.61 


0.006 






432 


2.61 


-50.52 


1, 


.27 


2.14 


0.77 


-44.29 


1, 


.31 


2.22 


2.21 


0.007 


BO 


0.50 


104 


2.25 


— 17.72 


i. 


.oy 


1.92 


0.15 


— 15.74 


i. 


.Do 


2.01 


2.00 


0.006 






225 


3.39 


-24.32 


1. 


.oy 


2.50 


0.71 


-20.89 


-1 

1, 


.Do 


2.67 


2.66 


0.010 






473 


5.23 


-34.83 


-1 

1. 


.oy 


3.22 


2.88 


-29.20 


-1 

1, 


.Do 


3.65 


3.12 


0.093 




0.70 


225 


3.39 


-34.05 


\ 


.27 


1.60 


0.39 


-29.25 


I 


.31 


2.08 


1.51 


0.334 






473 


5.23 


-48.76 


1, 


.27 


2.22 


7.38 


-40.88 


1, 


.31 


2.62 


1.57 


1.039 


B2 


0.50 


58 


4.79 


—6.47 


1 

1. 


.oy 


2.14 


2.07 


—5.44 


1 

1, 


.Do 


2.53 


2.10 


0.271 






121 


7.46 


-9.54 


1 

1. 


.oy 


2.64 


8.02 


-7.88 


1 

1. 


.Do 


3.07 


2.16 


0.969 






247 


11.70 


-14.44 


1 
1 


.oy 


3.16 


31.18 


-11.80 


1 

1, 


.Do 


3.65 


2.20 


3.542 




0.70 


121 


7.46 


-13.35 


1 
1 


01 
.A 1 


1.70 


16.22 


-11.04 


1 

1 , 


.ol 


2.00 


1.37 


7.580 






247 


11.70 


-20.22 


\ 


.27 


2.26 


63.36 


-16.52 


l_ 


.31 


2.60 


1.37 


29.788 


B5 


0.50 


62 


15.13 


—3.76 


-1 

1. 


.oy 


2.94 


68.00 


—3.06 


-1 

1, 


.Do 


3.39 


1.96 


12.335 






125 


23.89 


—5.84 


1 


.oy 


3.54 


266.14 


—4.72 


-1 

1, 


.Do 


4.07 


2.00 


42.823 






253 


37.80 


-9.13 


-1 
i 


.oy 


4.24 


1049.25 


-7.35 


i, 


.Do 


4.85 


2.00 


170.223 




0.70 


62 


15.13 


-5.27 


-1 
i 


ol 
.Z i 


2.10 


119.78 


-4.28 


i, 


Q1 


2.40 


1.35 


74.824 






125 


23.89 


-8.17 


1 
1 


07 
.Z / 


2.64 


474.30 


-6.60 


1 

i , 


1 

.oi 


3.00 


1.35 


296.062 






253 


37.80 


-12.79 


\ 


.27 


3.22 


1885.03 


-10.30 


I 


.31 


3.66 


1.35 


1175.860 


B9 


0.50 


63 


26.30 


—3.08 


1 
1 


.oy 


3.36 


290.95 


—2.49 


1 

1 , 


.Do 


3.85 


1.90 


61.286 






126 


41.63 


—4.83 


1 
1 


.oy 


4.04 


1146.13 


—3.88 


1 

1 , 


.Do 


4.61 


1.92 


224.1299 






254 


65.96 


-7.60 


1 
1 


.oy 


4.80 


4540.98 


-6.10 


1 

1 , 


.Do 


5.49 


1.92 


893.279 




0.70 


63 


26.30 


-4.31 


1 
1 


.z / 


2.48 


479.826 


-3.48 


1 

1 , 


.ol 


2.84 


1.33 


4867.086 






126 


41.63 


-6.76 


-1 
1 


.z / 


3.06 


1905.91 


-5.44 


-1 

1, 


.oi 


3.48 


1.35 


1177.084 


(X Ara 


0.50 


57 


3.97 


—4.29 


1 


.59 


3.026 


0.67 


—3.64 


1, 


.68 


3.21 


3.20 


0.004 






119 


6.15 


—6.22 


1 


.59 


4.16 


2.68 


—5.18 


1, 


.68 


4.45 


4.26 


0.026 






245 


9.63 


—9.33 


1 


.59 


4.92 


10.61 


—7.66 


1, 


.68 


5.81 


4.30 


0.202 




0.70 


57 


3.97 


-6.00 


1 


.27 


2.64 


1.96 


-5.10 


1, 


.31 


2.82 


2.31 


0.048 






119 


6.15 


-8.71 


1 


.27 


3.08 


7.63 


-7.26 


1, 


.31 


3.64 


2.37 


0.209 






245 


9.63 


-13.06 


1 


.27 


3.60 


29.84 


-10.72 


1, 


.31 


4.20 


2.45 


0.765 


a Eri 


0.50 


60 


7.97 


-2.77 


1 


.59 


4.44 


5.88 


-2.28 


1, 


.68 


5.21 


3.98 


0.126 






124 


12.52 


-4.20 


1 


.59 


5.06 


23.22 


-3.43 


1, 


.68 


5.93 


4.04 


0.535 






251 


19.75 


-6.49 


1 


.59 


5.84 


91.96 


-5.26 


1, 


.68 


6.79 


4.12 


2.044 




0.70 


60 


7.97 


-3.87 


1 


.27 


3.20 


16.24 


-3.19 


1, 


.31 


3.74 


2.17 


0.604 






124 


12.52 


-5.88 


1 


.27 


3.78 


63.68 


-4.80 


1, 


.31 


4.36 


2.25 


2.122 






251 


19.75 


-9.08 


1 


.27 


4.44 


251.31 


-7.36 


1, 


.31 


5.12 


2.33 


7.462 


a Col 


0.50 


62 


17.10 


-1.50 


1 


.59 


6.36 


54.50 


-1.22 


1, 


.68 


7.45 


5.02 


0.807 






126 


27.02 


-2.33 


1 


.59 


7.32 


216.70 


-1.88 


1, 


.68 


8.51 


5.10 


3.104 






253 


42.76 


-3.66 


1 


.59 


8.46 


863.29 


-2.94 


1, 


.68 


9.83 


5.18 


11.809 




0.70 


62 


17.10 


-2.10 


1 


.27 


4.80 


158.30 


-1.70 


1, 


.31 


5.58 


3.05 


2.243 






126 


27.02 


-3.27 


1 


.27 


5.62 


627.78 


-2.64 


1, 


.31 


6.50 


3.11 


8.402 






253 


42.76 


-5.12 


1 


.27 


6.56 


2497.14 


-4.12 


1, 


.31 


7.58 


3.17 


31.504 



Stellar models refer to those listed in Table [T] In these models, dipole b = 3 and wind /S^^j = 1. In all these models, To is 
of order 1 to 10 X lO^K, Us is in the range from 50 to 200 kms, t^^fi^^ ~ 1 to 40 mins. The values of tcoo; ^'''e of the order 
10 hours to a few days and depend on the properties of the stars considered. In the case of the last three stellar models, 
they are of the order of a few days to several months. 



the values of U shock are approximately 30% of the values of the wind flow speed Wn.w- An important feature of the results is 
that for given values of the rotation speed and surface magnetic field strengths, when the wind speed is faster the values of 
X2 are larger. This is because wind flows that are faster lead to greater heating of protodiscs causing earlier onset of MRI. 
This is apparent when we compare the protodisc properties for the theoretical model B2 star and the models corresponding 
to a Ara given in table [S] 

In stellar models of earlier types such as 06.5 and BO shown in table [S] we find that MRI occurs almost instantaneously 
so that the protodisc phase of the disc region ends almost as soon as it starts. In the case of later type stars such as B5 or B9 
shown in the same table, the time taken for MRI to occur can be several months or years depending on the field strength. For 
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Table 4. Values of the equatorial density, height 
and time taken when a point X becomes the 
outer edge in a protodisc of a uniformly rotating 
stellar model with the same parameters as a Ara 
with 5* = 0.7 and Bm,o = 245G. 



X 


PX.edge 




fix, edge/ R 


,edge 




(gm cm' 






(days) 


1.288 


l.lOE - 


12 


0.156 


29.836 


1.408 


2.03E - 


13 


0.164 


8.691 


1.608 


4.21E - 


14 


0.185 


3.502 


1.808 


1.34B - 


14 


0.208 


1.978 


2.008 


5.24E - 


15 


0.230 


1.278 


2.208 


2.32E - 


15 


0.251 


0.889 


2.408 


1.13E - 


15 


0.270 


0.648 


2.608 


5.87E - 


16 


0.288 


0.489 


2.808 


3.23E - 


16 


0.302 


0.378 


3.008 


1.86E - 


16 


0.314 


0.297 


3.208 


l.UE - 


16 


0.322 


0.237 


3.408 


6.90E - 


17 


0.325 


0.190 



A necessary condition for the onset of MRI is 
not satisfied when the rotation is uniform. 



Table 5. Values of the equatorial density, height and time taken when a point X 
becomes the outer edge or when MRI occurs at X in a protodisc of a differentially 



rotating 


stellar model with the same 


parameters as 


a Ara with 5*. eg = 0.7, 


e = 0.2 


and Bm 


,o = 245G. 














X 


PX,edge 




PX,MRI 




hx,edge/ R 


hx.MRI /R- ix,edge * 


X,MRI 




(gm cm" 




(gm cm' 






(days) 


(days) 


1.330 






2.62E - 


13 




0.149 


8.384 


1.410 






1.3AE - 


13 




0.158 


5.771 


1.610 






3.56E - 


14 




0.183 


2.973 


1.810 






1.27E - 


14 




0.207 


1.880 


2.010 






5A2E - 


15 




0.232 


1.328 


2.210 






2.61E - 


15 




0.256 


1.005 


2.410 






1.38E - 


15 




0.280 


0.797 


2.450 






1.22E - 


15 




0.284 


0.765 


2.510 


imE - 


15 






0.290 


0.704 




2.610 


7.35B - 


16 






0.300 


0.617 




2.810 


4.08B - 


16 






0.318 


0.481 




3.010 


2.37E - 


16 






0.333 


0.383 




3.210 


1A3E - 


16 






0.346 


0.309 




3.410 


8.92E - 


17 






0.355 


0.252 




3.610 


5.72E - 


17 






0.361 


0.207 




3.810 


3.77E - 


17 






0.362 


0.171 




4.010 


2.5iE - 


17 






0.357 


0.141 




4.190 


1.81E - 


17 






0.347 


0.118 





most of this time, the radial extents of such protodiscs will be small because the rate at which they shrink is faster initially 
and it slows down drastically as Xedge approaches Xi„„er. Note that there is a difference in the time taken for the onset 
of MRI between the protodiscs of these models and the model corresponding to the B7 star a Col, where MRI occurs in a 
shorter time. The reason for this is the fact that the wind speed for the a Col model is much faster, which causes greater 
shock heating of the protodisc. This is also the reason why there is a corresponding difference between the protodiscs of the 
B2 model and the model corresponding to a Ara in table [S] Protodiscs of stellar models with the properties of a Ara and a 
Eri do not show any significant differences. 

Tables |4] shows critical values of properties of a protodisc when a point X on the equatorial plane becomes the outer 
edge at X for a uniformly rotating stellar model having the properties of a Ara with 5* — 0.7. Tables [5] and [6] show the 
corresponding values when X becomes the outer edge or when MRI occurs at X for differentially rotating models with 
5*, eg ~ 0.7 and Si,, eg ~ 0.5, respectively, e = 0.2 and having the same basic stellar properties. Suppose that we consider the 
same basic model with different rotation rates and different magnetic field strengths. Then, when the rotation rate is fixed and 
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Table 6. Values of the equatorial density, height and time taken when a point X 
becomes the outer edge or when MRI occurs at X in a protodisc of a differentially 



rotating 


stellar model with the same parameters as a Ara with 5* en = 0.5, e = 0.2 


and Bm 


o = iiyCj. 








X 


PX,edge 


PX,MRI 


hx,edge/R fl-X,MRl/R 


tx,edge tx,MRI 




(gm cm ^) 


(gm cm ^) 




(days) (days) 


1.70 




5.34B - 15 


0.163 


0.577 


2.00 




1.37E ~ 15 


0.188 


0.322 


2.40 




3.40E - 16 


0.215 


0.186 


2.80 




\.\2E - 16 


0.235 


0.122 


3.20 




5.27E - 17 


0.242 


0.091 


3.60 




2.00E - 17 


0.234 


0.059 


4.00 




9.95E - 18 


0.201 


0.039 


4.20 




7.24E - 18 


0.167 


0.029 


4.24 




6.81B - 18 


0.159 


0.026 


4.26 




6.61E - 18 


0.154 


0.025 


4.28 


6.39E - 18 




0.148 


0.024 


4.38 


5.35E - 18 




0.099 


0.015 


4.44 


4.83E - 18 




0.043 


0.006 



Table 7. Variation of properties of protodiscs with changes in rotation speed and magnetic field strength for 
a stellar model with the parameters of a Ara, dipole 6 = 3 and ji^el = 1 for uniform and differential rotation. 



Uniform Rotation with <S* = Si,,eq Differential Rotation with e = 0.2 



S*.eq Bm.o ^^Zim -^(^,o -^-inner -^l ^shrink ^(f>,o ■^-inner -^l -^2 ^2 

(G) (G) (days) (G) (days) 



0.40 


245 


9.63 


-7.46 


1.84 


6.05 


6.19 


-6.13 


1.97 


6.62 


5.79 


0.089 


0.50 


245 


9.63 


-9.33 


1.59 


4.92 


10.61 


-7.66 


1.68 


5.81 


4.30 


0.202 


0.60 


245 


9.63 


-11.20 


1.41 


4.16 


17.89 


-9.19 


1.47 


4.88 


3.27 


0.378 


0.70 


245 


9.63 


-13.06 


1.27 


3.60 


29.84 


-10.72 


1.31 


4.20 


2.45 


0.765 


0.80 


245 


9.63 


-14.93 


1.16 


3.17 


47.54 


-12.25 


1.19 


3.68 


1.55 


3.690 


0.90 


245 


9.63 


-16.80 


1.07 


2.82 


65.38 


-13.79 


1.08 


3.27 


1.08 


70.030 


0.75 


57 


3.97 


-6.43 


1.21 


2.40 


2.51 


-5.47 


1.24 


2.73 


1.84 


0.114 


0.75 


119 


6.15 


-9.34 


1.21 


2.86 


9.74 


-7.77 


1.24 


3.35 


1.94 


0.387 


0.75 


245 


9.63 


-14.00 


1.21 


3.38 


38.05 


-11.49 


1.24 


3.93 


2.04 


1.287 


0.75 


499 


15.15 


-21.43 


1.21 


3.98 


149.80 


-17.43 


1.24 


4.59 


2.12 


4.453 



the magnetic field strength is increased, the value of X2 increases slightly. The results for stellar models with the properties of 
a Ara are shown in table[7]for uniform rotation and for differential rotation with e — 0.2. Also, when the stellar rotation rate 
is increased while the magnetic field strength is fixed, the value of X2 decreases. We note that the values of X2 for models with 
differential rotation are somewhat larger than those for the corresponding models with uniform rotation. Our results show 
that when stars have solar-type differential rotation, there are broad ranges of stellar magnetic field strengths and rotation 
rates such that 5*, eg 0.9 for which MRI sets in over a significant region from Xirmer to X2 before the protodiscs shrink 
radially to that region. However, when a star rotates with near critical rotation speed such that 5* > 0.9, we find that the 
time for the onset of MRI is longer than the time for the protodisc region to shrink to a thin ring, unless some other MHD 
instability occurs. 



8 EVOLUTION OF DISC MATERIAL AFTER ONSET OF MRI 

A detailed study of the evolution of disc material after the onset of MRI is beyond the goals of the present paper. However, 
we consider several factors that will play a role in the evolution from the protodisc phase to a quasi-steady disc. Because the 
electrical conductivity is extremely high in a protodisc region, the magnetic diffusivity is negligible and the Ohmic diffusion 
time is extremely long ( Maheswararj[20o3 ) . Thus, before the onset of MRI, the magnetic field within a protodisc region remains 



frozen in the material as it rotates. At the outer edge, the sum of the centrifugal force and thermal pressure exceeds the sum 
of gravity and the magnetic force, so that disc material will flow outwards and draw out magnetic field lines with it. We focus 
on stellar models where magnetic fields are sufficiently strong to channel wind fiow into disc regions. In the case of B type 
stars, field strengths of order IG to lOG would be appropriate. We wish to understand how the onset of MRI in protodiscs 



14 M. Maheswaran and J. P. Cassinelli 



of such stars helps in the formation of quasi-steady discs that survive over a period of several decades. Turbulence arising 
from MRI may persist for a long time and will cause changes in the magnetic field and the magnetic diffusivity of the disc 
material will be affected. The occurrence of MHD dynamo action is possible. A study of the evolution of the magnetic field 
will be complicated. In a quasi-steady state we expect that the structure and strength of the field will allow a slow outflow of 
the disc material. 

After the onset of MRI, viscosity becomes important and the Shakura-Sunyaev viscosity coefficient is given by q « 
B^/{Anpc^) subject to the condition that a ^ 1. The azimuthal viscous force is negative and will tend to decelerate the disc 
rotation speed. Also, in general, dB^/dz will no longer be zero and the azimuthal magnetic force may play a r ole in the 
rotat i on of the disc material. The importance o f viscosity has been invoked in the modelling of decretion discs (e.g.. iLee et a] 
I991I : lOkazakil I2OOII ) and accretion discs (e.g., IPringlel Il98ll : iBalbusI [iooil ) . In decretion discs, material is fed into the disc 



region from the stellar equator and drifts outwards due to the effect of viscosity. In accretion discs, material in Keplerian 
discs fall towards the central star due to viscous loss of angular momentum. These studies consider thin discs and set up 
equations for physical quantities averaged over the heights of discs. We have not used such an approach in our models of 
protodiscs. However, because the heights of protodiscs are relatively small at the time that MR I occurs, thin disc models 
may be a convenient alternative in a study of the evolution of disc material after the onset of MRI. lHawlev et a ] (|l995l . ll996l ) 
have used a shearing box model for thin accretion discs and found that there will be an outward flux of angular momentum. 
We are interested in scenarios in which the angular velocity distribution of the disc material readjusts so that a quasi-steady 
disc may be formed, possibly with inflow and slow outflow. A model conforming to such a description has been developed by 
Brown et a] 1 20081 ) where wind material is fed into a disc region by a magnetic field and fiows slowly outwards with quasi- 



Keplerian rotation speed. Another possible quasi-steady model may be a magnetic rotator decretion disc model in which a 
protodisc that suffers MRI will provide material with super-Ke plerian ro t ation spee d to a decretion disc region. Such a model 
will be more realistic than the decretion disc model studied bv lLee et all (1991) and Okazaki ( 200 ll ) because the nonmagnetic 
models of Lee et al and Okazaki require the extremely special condition that Keplerian rotation must occur at the stellar 
equator. 

In the case of a quasi-steady disc, continuity of mass and magnetic fiux require that diu? puT^) / dzo = 0, divS = and 
an appropriate MHD induction equation be satisfied. The radial component of the momentum equation for axially symmetric 
quasi-steady discs with u^, = Q can be written as 



azu 
where 



_1_ 

47r 



"-0 



GM 
^2 



'dB^ 



1 dp , 
p dw 



dz 



dB, 
dvj 



Bj, 



_d_ 

dzu 



vjB^) 



(34) 



(35) 



is the radial component of the magnetic force. The thermal pressure term in equation (|34p may be negligible in a quasi-steady 
state where the disc material has cooled down and a radiative force may be included when appropriate. We are interested in 
quasi-steady models in which the radial speed is small an d the rotation speed is approximately equal to the magnetically 
modified Keplerian speed = [GM{1 — (^)/zu]^''^ defined in lMaheswara3 l|2003l ). where = vj^T-^j (GMp) is the ratio of the 
radial component of the magnetic force to gravity and is small. The azimuthal component of momentum satisfies 
1 d 



■cu azu p 



(36) 



'i-Kpvu'^ dvj 

The last term in this equation represents the a-viscosity force and will be negative. When u^^ is positive, the excess angular 
momentum brought in by the wind materia l is removed by the outward drift of the disc material. We expect that u-^ 
will be small as shown by Brown et al ( 20081 ) . This is because the continuity of mass in a quasi-steady state requires that 
PwVn.w = PpUn.D ac ross the shock boundary. If the disc region is isothermal in a quasi-steady state, we have pn ~ pw^n^w /^s 
I Cassinefli et alll2002l ) so that we obta n u-^^j « v„,w{c^/Vn^w)- Because Cs <C v„,w for the types of stars that we consider, we 
find that u^^ <^ Vn,w- The z-component of the momentum equation remains the same as (|24|) . In a detailed study we will 
have to include appropriate jump conditions at the shock boundary of the disc. 

At the inner end of a protodisc, if the azimuthal viscous force is sufficiently strong to decelerate the rotation speed, disc 
material will tend to fiow back towards the star as in the case of an accretion disc. However, this can compress the field 
lines of a dipole type magnetic field causing an increase in the magnetic pressure that will oppose the backward fiow. Also, 
continuing infiow of wind material will supply angular momentum and it may be possible for a quasi-steady state to exist 
for some length of time despite any backward fiow that may occur at the inner boundary of the disc region. Obviously, a 
detailed study will be required to determine the nature of the dynamics in this region. Although there may be other types of 
MHD instability that could affect the angular velocity distribution in a protodisc, we focus on MRI because we are able to 
determine the conditions that must be satisfied by the central star for its occurrence in a protodisc. An important aspect of 
our study of protodiscs of differentially rotating stars is that scenarios may exist in which MRI can intervene to redistribute 
angular momentum of the disc region and, possibly, permit the formation of quasi-steady disc regions. 
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9 DISCUSSION 

In this paper, we set up equations and jump conditions to study the formation, structure and evolution of protodiscs around 
rotating magnetic stars with winds. We obtain solutions for a variety of stellar models with different rotation rates and 
magnetic field strengths. When the magnetic field of a rotating hot star is strong enough to channel the flow of wind material 
towards the equatorial plane and field lines from the stellar surface cross the equatorial plane and return to the star, a protodisc 
with super-Keplerian rotation and no radial motion will be formed in an equatorial region. Initially, the sum of the magnetic 
force and gravity in this region is sufficient to prevent any radial outflow that may be caused by the centrifugal force and 
thermal pressure. As wind material continues to flow into the protodisc, its equatorial density and height increase with time. 
The resulting increase in centrifugal force relative to the magnetic force causes the material at the outer edge of the protodisc 
to move outwards so that the radial extent of the protodisc region decreases. Shock heating causes protodisc temperatures rise 
to values of order 1 to 10 x lO^K. Because the radiative cooling time is longer than the heating up time, protodiscs retain their 
high temperatures. There is no meridional flow of material within a protodisc and the continuing increase in density causes 
the thermal pressure in the protodisc to increase in relation to the magnetic pressure. When the rotation rate of the central 
star satisfies the condition 5*, eg ^ 0.9 and the angular velocity of the protodisc region decreases radially outwards, MRI will 
occur before the radial extent of the protodisc can decrease significantly. In the case of most stars with near critical rotation 
rates such that <S* > 0.9, the time taken for the onset of MRI is longer than the time by the protodisc to shrink to a thin ring. 
The evolution of protodiscs is strongly influenced by the inflow of wind material, centrifugal outflow and the occurrence of 
MRI. The protodisc phase of a disc region terminates either when it shrinks completely or when some instability such as MRI 
sets in. We derive conditions for centrifugal outflow at the outer edge of a protodisc and for the onset of MRI in a protodisc. 
We perform numerical computations to study the structure and properties of protodiscs of several different stellar models and 
to determine how the protodisc phase terminates. Our results show that there are situations involving different fleld strength 
domains that lead to the formation of protodiscs that suffer MRI soon after their formation or shrink to a thin ring when no 
instability intervenes. 

Equation (|32[) gives the angular velocity of a protodisc in terms of that of the central star. When the central star has 
uniform rotation, the protodisc has approximately uniform rotation and if no instability intervenes, its radial extent continues 
to decrease. A protodisc region may take a long time to shrink completely and during most of this time it would appear 
like a ring. However, it is not clear whether steady uniform rotation is possible in the surface regions of rapidly rotating 
early type stars because large scale meridional circulations must occur in the radiative envelopes of such stars. Meridional 
circulation will transport angular momentum and disrupt the angular velocity distribution to cause differential rotation. 
Studies by Maheswaran (1968) showed that steady rotation may be possible in the interior of the radiative envelope with 
angular velocity constant along field lines and having about 3 percent variation from uniform rotation. However, because of 
the sensitivity of meridional circulation near the stellar surface to perturbations in angular velocity, steady uniform rotation 
may not possible in the surface region. Other effo rts to develop mode ls of steady uniform or differential rotation in early type 



may not possible m tne surtace region, utner etto rts to aevelop moae ls ot steagy unitorm or aitterential rotation m early type 
stars with magnetic fields have been discussed by Strittmatteij ( 1970l ) and TassouJ ( 2000l l . We are not aware of any conclusive 



theoretical or observational evidence that early type stars with fast rotation and magnetic fields will have steady unform 
rotation along the surface. 

When the stellar surface has differential rotation with angular velocity increasing towards the equator, equation (|32p 
shows the angular velocity of the protodisc decreases radially outwards. In this case, MRI occurs in the protodisc within a 
few hours or days. We know that su ch rotation occ urs in the Sun, where th e differ ential rotation is maintained by meridional 
circulation and convection (see e.g., Tassoull 2000l l. Recently, Maeder et all ( 2008h have suggested that convection zones will 



be present near the surface of rotating massive stars. This situation is the same as in the surface regions of the Sun where 
convective motions and meridional circulation are present. Thus, we suggest that it is reasonable to expect that there are 
massive stars with solar-type differential rotation. 

Our study of protodiscs was motivated in large part by t he cla im that discs will not form around magnetic rotator 
stars because of ce ntrifugal breakout found by lud Doula et al liooi) in numerical simulations of uniformly rotating stars. 



Maheswaran (|2003l ) had previously pointed out that centrifugal outflow of protodisc material due to increased density could 



diffuse into Keplerian, especially when viscosity was present. In studying the evolutionary properties of protodiscs, we have 
focussed on centrifugal outflow and the onset of MRI. It is possible that there are other MHD instabilities that may also play 
a role in the evolution of protodiscs to quasi-steady discs. We have been able to show that under certain highly plausible 
conditions at the stellar surface, MRI can occur in protodiscs before they shrink significantly through centrifugal outflow. 
Protodiscs in which MRI occurs in an inner part before the outer edge shrinks to that part are of interest because viscosity 
arising from the onset of MRI may be able to redistribute the angular velocity of the disc so that material may diffuse into a 
quasi-steady disc region with slow outflow. Although, constructing detailed models of quasi-steady discs is beyond the scope 
our present work, we simply point out that there are such models that are plausible. In the first model, the disc magnetic 
field remains connected to the stellar field and a quasi-stead y disc forms with infiow from the wind zone and slow radially 
outward flow with fleld lines stretching out with the flow (e.g., Brown et al 20081 ). Another possible quasi-steady model is one 
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in wh ich a protodisc evolves after the onset of MRI into a magnetic rotator decretion disc (analogous to the model of lLee et al 



1991 ) with inflow of sup er-Keplerian wind material. Therefore, the centrifugal breakouts found in the numerical simulations 



of lud Doula et all l|2006f ) in uniformly rotating models does not mean that quasi-steady discs cannot form around differentially 
rotating magnetic rotator stars. 

We have not addressed the question of how protodiscs may be detected observationally. Because the duration of the 
protodisc phase in most cases is rather short, being of the order of a few hours to many days, special observational efforts will 
be required during disc formation phases to infer their presence. In some special cases, protodiscs can survive for a few years 
in the form of rings. 
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